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^ ! Abstract 
(N 

i-H 

|Q ■ The polarization of inclusive T(nS) at the Fermilab Tevatron is calculated 

' within the nonrelativistic QCD factorization framework. We use a recent 

a: 

O 1 determination of the NRQCD matrix elements from fitting the CDF data on 

' bottomonium production from Run IB of the Tevatron. The result for the 



polarization of T(15) integrated over the transverse momentum bin 8 < pt < 



5— i 

■ 20 GeV is consistent with a recent measurement by the CDF Collaboration. 

The transverse polarization of T(1S) is predicted to increase steadily for px 
greater than about 10 GeV. The T(2S) and T(35) are predicted to have 
significantly larger transverse polarizations than T(IS'). 

PACS numbers: 13.85.t, 13.85.Ni, 14.40. Gx 
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The production of heavy quarkonium provides an ideal testing ground for our under- 
standing of the production mechanisms for heavy quarks and the nonperturbative QCD 
effects that bind the heavy quark- ant iquark pair into quarkonium. The purely leptonic de- 
cays of the J PC = 1 quarkonium states allow very clean measurements of their cross 
sections. One observable that is particularly sensitive to both the heavy-quark production 
mechanisms and the QCD binding effects is the polarization of the J PC = 1 states, which 
can be measured from the angular distribution of the leptons from their decays. 

The NRQCD factorization approach provides a systematic framework for calculating the 
inclusive production rates of heavy quarkonium 0. It is based on organizing the inclusive 
production rate into a double expansion in powers of a s (mq), the running coupling constant 
at the scale of the heavy-quark mass, and of vq, the typical relative velocity of the heavy 
quark in the quarkonium rest frame. It should therefore be more reliable for bottomonium 
than for charmonium, because both of the expansion parameters a s (mg) and vq are smaller. 
In the NRQCD factorization approach, the larger-than-expected cross sections observed at 
hadron colliders are explained by introducing phenomenological parameters that describe 
the probabilities for the formation of the quarkonium state from color-octet heavy quark- 
antiquark pairs Once those parameters are determined, the cross sections for polarized 
states are predicted without any additional parameters. 

The most dramatic prediction that has emerged from the NRQCD framework is that 
in hadron colliders, J PC = 1 quarkonium states should become increasingly transversely 
polarized as the transverse momentum increases f3|]4[]. Quantitative predictions of the 
polarization of the charmonium states J/ijj and ip{2S) at the Fermilab Tevatron indicate that 
there should be little polarization for pt around 5 GeV, but that the transverse polarization 
should increase steadily at larger px |§H^I- The first measurements of the polarization of 
J/ijj and ip(2S) by the CDF collaboration || have shown no evidence for this predicted 
increase. However, the discrepancies with the theoretical predictions are significant only in 
the highest pt bin, so a definitive conclusion must await the higher statistics measurements 
that will be possible in Run II of the Tevatron. 
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The CDF Collaboration has recently measured the polarization of inclusive T(IS') in 
Run IB of the Tevatron ||. The results for the px bins from 2 to 20 GeV and from 8 
to 20 GeV are both consistent with no polarization. Since the cross section falls rapidly 
with pt, this indicates that there is little if any polarization for pt below about 10 GeV. 
To determine whether this result is compatible with the NRQCD prediction, we need a 
quantitative calculation of the polarization for inclusive T(1S) mesons. 

The theoretical ingredients needed to calculate the polarization of direct T(1S) mesons 
have been available for several years f6|jl(l]. They were used by Beneke and Kramer and by 



Leibovich to predict the polarization of prompt ip(2S) at the Tevatron |5|||. The calculation 
of the polarization of inclusive T(nS) is complicated by the fact that the inclusive signal 
includes T(1S) mesons that come from decays of higher bottomonium states. Decays of 
XbO-P), T(2S), and x&(2-P) account for about 27%, 11%, and 11% of the inclusive T(1S) 
signal, respectively [[11]]. The missing ingredients in the calculation of the polarization of 
inclusive T(IS') were the cross sections for polarized Xbj- The necessary parton cross sections 
were recently calculated by Kniehl and Lee and used to predict the polarization of 



prompt J /if) at the Tevatron |7]]. 

In this paper, we present quantitative predictions for the polarization of inclusive Y(IS'), 
T(2S), and Y(3>S) at the Tevatron using the NRQCD factorization formalism. We use 
NRQCD matrix elements for direct bottomonium production which were recently determined 
by Braaten, Fleming, and Leibovich from an analysis of data from Run IB at the Tevatron 
|13fl . Our result for the polarization of T(IS') is consistent with the recent measurement 



by the CDF Collaboration in the pt bin from 8 to 20 GeV. We predict that the transverse 
polarization of T(IS') should increase steadily for pt greater than about 10 GeV, and that 
the T(2S) and T(3S) should be even more strongly transversely polarized. 

Since the T(nS) is a spin-1 particle, the projection A of its spin along any quantization 
axis should be —1, 0, or +1. The polarization of T(nS) can be measured from the angular 
distribution of the leptons from its leptonic decays. The angular distribution of the positive 
lepton with respect to the Y(nS) momentum in the hadron center-of- momentum frame is 



proportional to 1 + a cos 2 8, which defines the polarization variable a. Taking the T(nS) 
momentum in this frame to be the spin quantization axis, we define the longitudinally 
polarized T(nS) to be the A = state and denote it by T^(nS'). The polarization variable 
a is related to the fraction £ of T(nS) mesons that are longitudinally polarized by a = 
(1 — 3£)/(l+£). The longitudinal polarization fraction £ = jo is the ratio of the inclusive 
cross section for Tl(uS) to the inclusive cross section for T(nS) summed over spins. 

The cross section a for inclusive T(nS) is the sum of the direct cross section and the 
direct cross sections for the higher bottomonium states T(mS) and Xb(mP) weighted by the 
inclusive branching fractions Bjj^r(nS) fc> r H — > T(nS) +X: 

a[T(nS)] inc = a[T{nS)} + £ B H ^ r(nS) a[H] . (1) 

H 

For T(3S), we consider only the direct channel, neglecting any possible feeddown from 
higher states, such as x&(3P). For T(2S), we take into account the direct channel and the 
feeddown from ^(2P) and T(3S'). For T(15), we take into account the direct channel and 
the feeddown from x 6 (lP), T(2S), x&(2P), and T(3S). The inclusive branching fractions 
B H _+~[(nS) are given in Table I of Ref. . 



The cross section a L for inclusive T L (nS) is the sum of the direct cross section for T L (nS) 
and the direct cross sections for the higher spin states T(mS)x and Xbj( m P)x weighted by 
B H ^ T {nS) an d by the conditional probability Ph x ^t l (uS) f° r H\ to decay into T L (nS) given 
that it decays into T(nS): 

a L [T(nS)} mc = a[T L (nS)]+J2 B H-+r(nS)PH x -+r L (nS)cr[Hx}, (2) 

The conditional probabilities are given in Table |. For Xbj(nP)\ — > Ti(nS), n = 1,2, they 
are given by simple angular-momentum factors for the radiative transition. For Xbj{2P)\ 
T L (1S) and T(mS)\ — > T L (nS), we must average over the various decay paths weighted 
by their branching fractions. The important steps in the decay paths are of 3 kinds. The 
observed hadronic transitions T(mS) — > T(nS) + tot preserve the spin A. For the radiative 
transitions Xbj(2P)\ Tl(L9) + 7 and T(mS)\ — > Xbj{ n P)\' + 7 , the probabilities for 



each spin state are given by simple angular-momentum factors |[L4|| . Taking the weighted 
average over the decay paths, we obtain the results in Table |. 

The polarization variable a for T(nS) has been expressed as a ratio of linear combina- 
tions of the direct cross sections for T(nS) and higher bottomonium states. The NRQCD 
factorization formula for the direct cross section for a bottomonium state H of momentum 
P and spin quantum number A has the schematic form 

da\pp^H x (P)+X] = ^2da\pp^bb n (P)+X}(0^), (3) 

n 

where the summation index n extends over all the color and angular momentum states of 
the bb pair. The bb cross section can be expressed as 

da\pp -»■ bb n (P) + X] = f i/p ® f j/p ®da[ij ^bb n (P)+X], (4) 

where fi/ p (x,fi) and fj/p(x,fi) are parton distribution functions (PDF's) and a sum over the 
partons i,j is implied. The parton cross sections da can be calculated using perturbative 
QCD. All dependence on the state H is contained within the nonperturbative matrix ele- 
ments (0^ x ^). In general, they are Lorentz tensors that depend on the momentum P and 
the polarization tensor of H\. The Lorentz indices are contracted with those of da to give 
a scalar cross section. The symmetries of NRQCD can be used to reduce the tensor matrix 
elements (0^ x ^) to scalar matrix elements (0%) that are independent of P and A. This 
reduces the variable a to a ratio of linear combinations of the NRQCD matrix elements. 

A nonperturbative analysis of NRQCD reveals how the various matrix elements scale 
with the typical relative velocity v of the heavy quarks. The most important matrix el- 
ements for the production of the S-wave states T(nS) and rjb(nS) can be reduced to one 
color-singlet parameter (Oi^ nS \ 3 Si)) : which scales like t> 3 , and three color-octet parame- 
ters (Oj {nS \ 3 S 1 )), (Oj (n5) (%)>, and (Oj {nS) { 3 P )), all of which scale like v 7 . The most 
important matrix elements for the production of the P-wave states Xbj{nP) and h{nP) 
can be reduced to a color-singlet parameter (Of b0 ( 3 Po)) an d a single color-octet parameter 
(Og b0 ( 3 Si)) , both of which scale like v 5 . At higher orders in v, so many new matrix elements 
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0.853 ± 0.009 




±1 
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±1 
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TABLE I. Conditional probabilities Ph x -*t L (nS) f° r the bottomonium spin state H\ to decay 
into a longitudinally polarized T(nS) given that it decays into T(nS). 
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enter that the predictive power of the NRQCD approach is lost. We therefore assume the 
matrix elements enumerated above are sufficient to describe the bottomonium cross sections. 
The first determination of the color-octet matrix elements for bottomonium production 



was a pioneering analysis by Cho and Leibovich [IS] of the data on bottomonium production 



from Run IA of the Tevatron [IB]. Due to the limited statistics, they had to use educated 



guesses for some of the matrix elements. An updated theoretical analysis based on the 
new CDF data from Run IB |§ has been made by Braaten, Fleming, and Leibovich Jl3| . 



Their color-singlet matrix elements are given in Table II of Ref. []13| . Those for T(nS) 
were determined phenomenologically from the leptonic decay rates of T(nS), while those for 
Xbj{nP) were estimated from potential models. Their color-octet matrix elements are given 
in Table V of Ref. [|13[j. They were determined by fitting the CDF data on the differential 
cross sections for T(IS'), T(2S), and T(3S) at px > 8 GeV and from the fractions of T(1S) 
from the decays of Xb(l-P) and x&(2-P) [fTTjl . 

The leading terms in the parton cross sections in (|^) depend on the region of the trans- 
verse momentum px of the bb pair. For px in the range 8 GeV < pt < 30 GeV, the leading 
terms are fusion contributions from the parton processes ij — > bb + k. For px much greater 
than 2mft, fragmentation contributions from parton processes such as ij — > g + k, followed by 
g — > bb, become important. For charmonium, fragmentation effects change the differential 
cross sections by less than 3% at px = 5 and less than 11% at px = 10 GeV. Since m& 
is 3 times larger than m c , we expect fragmentation effects to change the differential cross 
sections for bottomonium by less than 11% for px less than 30 GeV. We will therefore avoid 
the complications of fragmentation by restricting our predictions to pt < 30 GeV. 

For pt much smaller than 2m>b, parton processes such as ij — > bb + ggg.-- involving the 
multiple emission of soft gluons become important and it is necessary to resum their effects. 



In the analysis of Ref. ||13|| , this problem was avoided by using only the CDF data for px > 8 
GeV to fit the color-octet matrix elements. Since we will be using the matrix elements from 
that analysis, we will also restrict our predictions to px > 8 GeV. 

Having restricting our attention to the region 8 GeV < px < 30 GeV, it should be safe to 



use only the fusion cross sections for da in (|j). We include the parton processes ij — > bb + k, 
with i,j,k = g,q,q and q = u,d,s,c. We treat the c quark as a massless parton. The 
leading-order parton cross sections da are proportional to a^(/i). They are given in Refs. pfl 



and pll for all the relevant bb color and spin states with the exception of color-singlet Pj 



states, for which they are given in Ref. 



We follow the analysis of Ref. [ 13j as closely as possible. We use a common renor- 
malization and factorization scale \x for fi/ p , fj/p, and a s , taking its central value to 
be fix = (ml + Pt) 1 ^ 2 an d allowing it to vary within the range \[it to 2/it- We set 
rrib = 4.77 ±0.11 GeV. We consider two choices for the PDF's for comparison: CTEQ5L and 
MRST98LO [17]]. We evaluate a s (fi) from the one- loop formula with n/ = 5 using the bound- 
ary conditions a s (M z ) = 0.127 for CTEQ5L and a s (M z ) = 0.125 for MRST98LO. The cross 
section for direct T(nS) depends on the linear combination (O 8 ( 1 S , )) + r(O 8 ( 3 P )) /ml, with 
r varying from 4.6 at pt = 8 GeV to 3.2 at px = 30 GeV. In the cross section for direct 
Tl(uS), r varies from 7.6 at pt = 8 GeV to 3.6 at pt = 30 GeV. We consider two extreme 
cases: (Og^So)) = and (0$( 3 Po)) = 0. The NRQCD matrix elements and their statistical 
errors are given in Table II and V of Ref. |13[] . The color-octet matrix elements in Table V 
also have an additional upper and lower error associated with changing /x by a factor of 2. 
This allows the correlation between the errors in /i and the color-octet matrix elements to 
be taken into account. 

The polarization variable a can be expressed as a ratio of linear combinations of NRQCD 



matrix elements. The errors in the matrix elements from Ref. |13| give large uncertainties 
in the numerator and the denominator, but they tend to cancel in the ratio. As our central 
value for a, we take the average value from the 4 combinations corresponding to the two 
choices (O 8 ( 1 5 , )) = or (O 8 ( 3 P )) = and the two choices of PDF's. The deviations among 



1 There is a typographical error in Eq.(22) of Ref. ||i~2| . The color-octet matrix elements 
(0^°( 3 S{ 8) )) = (Cg*( 3 Si)) should be replaced by (2J + 1) (O^ ( 3 s{ 8) )). 
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the 4 combinations are included in our error. We compute the errors from the remaining 
parameters by varying each one independently for the choices MRST98LO and (Og( 3 Po)) — 
and using the central values of the other parameters. The parameters include the scale /x, 



the mass m&, 5 color-singlet matrix elements |]I3"||, 8 color-octet matrix elements ||13||, and 



18 exclusive branching fractions for quarkonium transitions [113. To compute the variation 
from /i, we change \i by a factor of 2 or | in the PDF's and in a s and simultaneously shift all 
the matrix elements by the second upper or lower error in Table V of Ref . fl3| . This takes 
into account the large correlation between the choice of fi and the values of the color-octet 
matrix elements. All the variations are added in quadrature to obtain the final error on a. 

Our results for the polarization variable a for T(IS'), T(2S), and T(3S) for the Tevatron 
at ^/~s = 2.0 TeV are shown as shaded bands in Figs. [1], @, and respectively. The curves in 
Figs. [l|-|3] are the central values of a for direct T(nS), for T(nS) from Xb(mP) + 7, and for 
T(nS) from T(mS) + 7nr. These channels together provide a complete decomposition of the 
inclusive rate. The fractions of the inclusive rate from each of these channels vary slowly 
with pt and add up to 1. For inclusive Y(IS'), the fractions have been measured by the CDF 
collaboration to be approximately 52% for direct T(IS') and 27%, 11%, and 11% for T(15') 
from x b (lP), T(2S"), and x fe (2P), respectively [TT|. For inclusive T(2S), the fractions can 



be estimated from the analysis of Ref. |13[] to be approximately 73% for direct T(2S') and 
25% and 2% for T(2S) from x b (2P) and T(3S), respectively. For T(35), we ignored any 
possible feeddown from higher states. 

The predictions for a for -y/i =1.8 TeV are essentially identical to the predictions for 
y/s = 2.0 TeV in Figs. [l|-|3|. The cross sections and a are both smaller by about 16%, 
but the change cancels in the ratio. Integrating over pt in the range 8 GeV < px < 20 GeV, 
we obtain a = 0.13 ± 0.18. This is in good agreement with the value measured by the CDF 
collaboration: a = 0.03 ± 0.28 @. 

In Figs. [1]-^], the curves for a for direct T(nS) and for T(nS) from T(mS) + irn increase 
steadily with p?. The curves for T(nS) from Xb(mP) + 7 increase at first, but then flatten 
out at a value of a just above 0.2. Thus the feeddown from \bijnP) tends to wash out the 
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polarization. The fraction of T(nS) from the Xb states is approximately 38% for n = 1, 
25% for n — 2, and probably smaller yet for n — 3. We therefore expect a to increase more 
rapidly with px for T(2S) and T(3S) than for T(IS'). A much larger transverse polarization 
for T(2S + 3S) than for T(IS') has also recently been observed in bottomonium production 
in p-Cu collisions at y/s = 38.8 GeV 

The predictions in Figs. [l]-[3] are based on NRQCD matrix elements extracted from CDF 
data in the range 8 GeV < pr < 20 GeV. The extrapolation of these results to lower values 
of pt should be considered unreliable. This is evident from the curves for direct T(IS') from 
T(2S) in Fig. [3] and direct T(2S) in Fig. [|. The dramatic changes in these curves at the 
smaller values of pt are artifacts of the fit of Ref. |13| having given negative central values 
for the color-octet matrix elements (O 8 ( 1 S , )) or (0 8 ( 3 -Po)) f° r T(2S). 

In Run II of the Tevatron, the much higher statistics will allow more accurate measure- 
ments of the polarization of T(1S) in several bins of pr- It may also allow measurements of 
the polarizations of T(25) and T(3S). Significant improvements in the theoretical predic- 
tions are also possible. The most important improvement is taking into account the effects 
of multiple soft-gluon emission at low px- This would allow more accurate determinations 
of the color-octet matrix elements, since the entire pt range of the CDF data from Run I 
could then be used in the fits. It is also important to include fragmentation effects, so that 
the predictions can be extrapolated with confidence to large px- 

In conclusion, we have calculated the polarization of the T(nS) states at the Tevatron 
using the NRQCD factorization approach. Our result is compatible with the recent CDF 
measurement for T(IS') in the pt bin from 8 to 20 GeV, which is consistent with no polar- 
ization. However, our results also indicate that a nonzero transverse polarization should be 
observable for all three T(nS) states in Run II of the Tevatron. 

J.L. thanks Sungwon Lee and Heuijin Lim for their useful advice. This work was sup- 
ported in part by the US Department of Energy Division of High Energy Physics under 
grant DE-FG02-91-ER40690 and by the KOSEF and the DFG through the German-Korean 
scientific exchange program DFG-446-KOR-113/137/0-1. 
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FIG. 1. Polarization variable a vs. pr at yfs = 2.0 TeV for inclusive T(IS') (shaded band). 
The curves are the central values for direct Y(IS'), Y(IS') from T(2S') + tttt, and T(15) from 
Xb(lP)+ 7 or X b(2P)+J. 
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FIG. 2. Polarization variable a vs. pr at y/s = 2.0 TeV for inclusive T(2S) (shaded band). 
The curves are the central values for direct T(2S), T(2S) from T(3S) + tttt, and Y(2S) from 
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FIG. 3. Polarization variable a vs. px at = 2.0 TeV for T(35) (shaded band). The curve 
is the central value. 
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